DOE/ER/40561-82-INT00-INT 



Extracting particle freeze-out phase-space densities and entropies 
from sources imaged in heavy- ion reactions 

David A. Brown 1 , Sergei Y. Panitkin 2 and George F. Bertsch 1 

O 

institute for Nuclear Theory, University of Washington, Box 351550, Seattle, WA 98195-1550 
Department of Physics, Kent State University, Kent, Ohio 44242 
■ (February 8, 2008) 

in 

<N 

0\ ; Abstract 
m 
o 

^ ■ The space-averaged phase-space density and entropy per particle are both 
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fundamental observables which can be extracted from the two-particle corre- 
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lation functions measured in heavy-ion collisions. Two techniques have been 
proposed to extract the densities from correlation data: either by using the 
radius parameters from Gaussian fits to meson correlations or by using source 



^ ■ imaging, which may be applied to any like pair correlation. We show that 
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the imaging and Gaussian fits give the same result in the case of meson inter- 
ferometry. We discuss the concept of an equivalent instantaneous source on 
which both techniques rely. We also discuss the phase-space occupancy and 
entropy per particle. Finally, we propose an improved formula for the phase- 
space occupancy that has a more controlled dependence on the uncertainty 
of the experimentally measured source functions. 
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I. INTRODUCTION 



The phase-space density of particles produced in an ultra-relativistic heavy-ion collision is 
a fundamental observable which is accessible, at least in part, via two-particle correlations. 
Measurements of this observable are interesting because they may either provide direct 
evidence for thermal phase-space distribution of particles at freeze-out or show evidence of 
deviations from such P],0]. Also, if one can measure the average phase-space density, one 
can begin to look for effect of an overpopulation of phase-space or more exotic phenomena 
such as pion lasers, superradiance, etc. M. A quantity closely related to the phase-space 
density is the entropy per particle, a key thermodynamic property of high density matter. 
Indeed, a phase transition might not take place adiabatically and could generate entropy. 

Several methods were suggested for the measuring the phase-space densities of the various 
particles in a heavy-ion collision. For the case of identical noninteracting pions, Bertsch Q 
proposed a method that uses both the radius parameters from conventional HBT analysis 
and the pion spectrum. He also suggested a physically intuitive definition of phase-space 
density: replace the single particle sources in the Koonin-Pratt formalism |5]|| with an 
instantaneous source at a mean freeze-out time. This effective source is just the effective 
phase-space density at freeze-out. This replacement preserves the particle numbers at large 
times (after freeze-out) even though, by construction, it differs from the true density during 
the freeze-out process. 

Another approach, based on the source imaging technique introduced by Brown and 
Danielewicz |7]|| , allows one to extract the space-averaged phase-space density and entropy 
per particle from any like-pair correlations. This approach uses both the source function 
imaged from a correlation measurement and the single particle spectrum. We will show that 
this approach is more versatile: it can be applied to any like-pair correlations and it reduces 
to the HBT result of when applied to pions. 

A third approach, introduced by Siemens and Kapusta ||, uses the ratio of total deuteron 
yield to total proton yield instead of the two-particle correlation function to estimate the 
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total phase-space occupancy of nucleons. This approach is complimentary to the source 
imaging approach in the sense that both give access to the proton phase-space occupancy. 
However, while the imaging method requires one to measure the correlation function and 
one particle spectrum for one specie of particle, the deuteron-proton ratio approach requires 
one to measure the spectra of two different species. Usually these spectra have different ac- 
ceptances and efficiencies of particle reconstruction and identification which may complicate 
experimental analysis. We will not discuss this approach in this letter. 

The outline of the paper is as follows. First, we use the substitution of Ref. @] to explain 
how the imaged sources can be used to extract the space-averaged phase-space density, 
(f(p))i as discussed in 0. We will show that these results for (f(p)) are a generalization of 
the result in Ref. |4| derived for identical mesons. Next, we will demonstrate that one does 
not change the source function by making this substitution. Within this framework, we will 
discuss the phase-space occupancy, (/), and the entropy per particle, S/N paTt . We will find 
an expression for the the phase-space occupancy that is an improvement over that in J7| as 
the formula here has a more controlled dependence on the uncertainty of the experimentally 
measured source functions. 



II. THE CORRELATION AND SOURCE FUNCTIONS 

We begin by defining a Lorentz-invariant single particle source, D(r,p), which gives the 
rate for creating on-shell (E = \ffp + m 2 ) particles (of all spin projections): 

, Ed 7 N 

In our notation, r = (t, r) is a four- vector and p is a three vector. With our choice of 
normalization, the single-particle source transforms as a Lorentz scalar. The single-particle 
source may be computed directly in an event generator such as RQMD [|1(| . 

We define the two-particle source function as the probability density for producing a 
particle pair separated by r' in their center of mass (CM) frame. Following we define 

this source function as the convolution of the single-particle sources: 



f dt' f d 4 RD(R + r/2,p 1 )D(R-r/2,p 2 ) 

n d^RD^RjS) (/ d 4 RD(R,p 2 )^j (2) 
= J dt' J d 4 RD(R + r/2,p 1 )D(R-r/2,p 2 ) 
Here we denote the normalized single-particle sources with a tilde, coordinates taken in the 
pair CM frame with primes, and coordinates taken in the the lab frame without primes. 
The average particle momentum is p = + p 2 ) and the relative particle momentum is 
q — \{pi — p 2 )- While all q pairs seem to enter in this definition (0), in practice only low 
relative momentum pairs contribute to correlations and hence either to the imaged sources 
or Gaussian fits to correlations. Typically, this is used to allow one to replace p\ and p 2 with 
the average pair momentum p (see Ref. 0). Following convention, we will not explicitly 
state the q dependence of the source function. 

The two-particle source function may be extracted from the measured two-particle cor- 
relation function by inverting the Koonin-Pratt [HH equation: 
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Here C$(q') is the measured correlation function, and $ -*, (r') is the pair relative wavefunction 
in the pair CM frame. We comment that the correlation function can be measured and the 
source function extracted in any frame, however relation between the two is simplest in the 
pair CM frame. 

In order to extract the source function, we first discretize Eq. (0) to obtain the matrix 
equation Cj = J^jKijSj- We then proceed as in [0,|J and find the set of source points, Sj, 
that minimize the x 2 - Here, \ 2 — J2i(Ci — J2j KijSj) 2 / A 2 Ci. The \ 2 minimizing source 
is Sj = J2i[{K 7 BK)~ l K 7 B]ji(Ci — 1) where K T is the transpose of the kernel matrix 
and B is the inverse covariance matrix of the data Bij = 5ij/A 2 Ci. The error on the 
source is the square-root of the diagonal elements of the covariance matrix of the source, 
A 2 ^ = (K T BK)~ l . Since this procedure works on any like-pair correlation, we can dispense 
with the correlation function and work directly with the source function. 



III. THE PHASE-SPACE DENSITY 



The phase-space density at time t is a Lorentz scalar; we write it as 

xl - -» 1 d&N 

Here we have defined a unit volume in phase-space as To = (2s + 1) / (2irhc) 3 and a differential 
phase space volume as T drdp. The (2s + 1) factor accounts for the spin of the particle of 
interest. Assuming that the particles propagate as free particles, the phase-space density at 
a specific time t\ can be written in terms of the phase-space density a different time t as 

f(t u f,p) = f(t ,f- v(t x - t ),p) 

1 ru (5) 



+— dtD(t,r-ti(t x -t),$) 
1 h/ J t 

This expression ignores processes that modify the particle's momentum (such as from a 
mean-field, Coulomb forces, etc.). Since we want the particles at position r, we must look 
backwards (or forwards) in time to where it was at its creation, e.g. r — v(t\ — t). 

A particle is said to have frozen out when it has undergone its last strong interaction and 
is now propagating to the detector. Since each particle can freeze-out at a different time, the 
concept of a "freeze-out phase-space density" is ambiguous. One option to deal with this 
is to average over the creation times of the particles. However then the question becomes 
how to perform this average and what the averaged density means. Ref. [§]] proposes a 
simple solution: replace the particle source D(r,p) with an effective source D eG (r,p) at an 
instantaneous freeze-out time, £/. In other words, 

£>( r> £) _> D cS (r,p) = 5(t f ~t)ET f cS (f,p) (6) 

The effective freeze-out phase-space density is related to the true single-particle source via: 

f cS (r,p) = ^ j\tD(t,f-v(t f -t),p) (7) 

Here, T > tf > —T and ±T is simply some large time after (before) which the source 
is turned off (on). It is clear from (|7|) that replacement in @ gives the right phase-space 



density at large times (t > T) as it must by Liouville's Theorem. Furthermore, this can be 
easily checked by examining Eq. (|5|) for large times t > T with / = for times earlier than 
—T. Finally, it is clear that this replacement does not alter the momentum space density 
that one might calculate from Eq. (|5|). 

Not only does this prescription give a reasonable way to define the freeze-out phase-space 
density, but it also is entirely consistent with our definition of the source function in Eq. (j^). 
We will illustrate this by simply inserting the effective phase-space density into the definition 
of the source (Eq. (g)) and performing the required algebra. Special care must be taken in 
performing the various time integrals as some are performed in the lab frame while others 
are performed in the pair CM frame. 

We now insert Eqs. (§,0) into the equation for the source function in Eq. (|2|) and perform 
one of the delta function integrals. We obtain 

S f (P) = J dt' 5{t) J d 3 R J dT J dr 

xD(T + r/2 1 R + r/2-vi(t f -T-r/2),p 1 ) (8) 
xD(T - r/2, R-f/2- v 2 {t f -T + r/2),p 2 ). 

Here Vi = Pi/Ei are the velocities of the individual particles. Also, T — (ti + t 2 )/2 and 
r — (t\ — t 2 ) are the average and relative time variables from the time integrals in the 
effective phase-space density. We can simplify Eq. (g) by introducing the average velocity of 
the pair: v = pj ^p 1 + m 2 so that = v + SHi. With this, we make the change of variables 
R — v(tf — T) — > R and remove much of the time dependence from the spatial arguments: 

Sp(?) = J dt' 6{t) J d 3 R J dT Jdr 

xD(T + r/2, R + (f + vr)/2 - 5v x {t f - T - r/2), ft) (9) 
xD(T- r/2, R-(r + vr)/2 -8v 2 (t f -T + r/2),p 2 ) 

Our next step is to do the if integral. This appears to be a straightforward delta function 
integral, however care must be taken due to the different reference frames involved. When we 
perform the t' integral, two things happen: / dt' 5{t) — > I/7 and everywhere in the integrand 
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t' — > — vrL Here r'u is the component of r' in the direction of the boost from the lab to the 
pair CM frame (this boost velocity is v). 

We make further progress by examining f+vr. First, define the four- vector s = (r, s) = 
(r, r+vr). If we write this in the pair CM frame, we have s = (y{r' + usj|), 7(sj| + vr'), r'j^j . 
So, the parallel component of s is 

T f 1 

S\\ = 7(s(i + VT 1 ) = 7(s(i + v( VSti)) = —s\i + VT (10) 

II II 7 II 7 

But note, from the definition of s and from the t' integral we already have 

sii = r\\ + vt = y(r\\ + vt') + vt = —r'u + vr. (11) 

ii 7 ii 

Therefore, if we identify sji = r'» and change integration variables from r to r', we find the 
result 

Sp(P) = Jdr' J d 3 R J dT 

xD(T + r/2, R + f/2- 5v x {t f - T - r/2),p x ) (12) 
xD(T- r/2, R-r/2-5v 2 (t f -T + r/2),p 2 ). 
In order to complete the connection between Eq. ( |12"D and the source in Eq. (fj) , we must 
justify dropping the Sv^s in Eq. fll2l) . Since At = (tf — T^ft/2) is on the order of the freeze- 
out duration, if SvtAt is smaller than the characteristic length scale of the single-particle 
source, we can drop the Sv^s. Writing 5vi in terms of the relative momentum q, we find 

6vi=v-Vi = ±hq-v-qv) + 0{q 2 ). (13) 
hi 

Thus, the 5vi term shifts the spatial argument of D in Eq. (|13D ~ Atq/yE in the direction 
parallel to v and ~ Atq/E perpendicular to it. For highly relativistic or massive pairs, this 
shift may be neglected. For low velocity or light pairs (such as tt pairs) the shift is important, 
especially if the freeze-out duration is large or the system size is small. 

A. Space- Averaged Phase-space Density from the Sources 

Inserting the particle source with instantaneous freeze-out into Eq. (^), taking pi,p 2 ~ p, 
taking the limit as r' —>■ 0, and performing the integrals over time, we find Eq. (17) of Ref. 0: 



Jd 3 rf 2 (r,p) 1 1 Ed 3 N 



(f(p)) = J — r = T^Sdr' -> 0) (14) 



For the sake of brevity, we drop both the "eff " tag on the density here and for the remainder 
of the paper. 

For zero impact parameter collisions, we may exploit the azimuthal symmetry and aver- 
age over the angle of the particle transverse momentum 8 PT and use 

f dB~. I d 3 r f 2 (r.v) 

(f(y,PT)) 



J d9p T 


[ d 3 rf 2 (f,p) 


J d0 PT 


1 d 3 rf(r,p) 



(15) 

1 1 d' 2 N „ 

S?(r> -> 0) 



r m 2ndy dp? Pt 

Note that, due to the azimuthal symmetry the density is a function of the particle rapidity 
y and the magnitude of the particle transverse momentum, px- 

As a side comment on using Eqs. ( |i"4"D and (|15|) , one can use either the r' — > point 

or 



from either a full three-dimensional reconstruction of the source function (such as in ||11|| 
from the angle-averaged source function. In practice, it is usually much easier to measure 
the angle-averaged two particle correlation function (and hence source function) then the 
full three-dimensional correlation because one can sum over angles to increase statistics. 

B. Gaussian Meson Sources 

We now show that Eqs. (|1^) and fll4] ) are direct generalizations of the results in M or 
for pions. The correlation function for identical non-interacting spin-0 bosons can be written 
in terms of a matrix of radius parameters ]T^] : 



C($) = l + \e- QiQ ' lB?]i *. (16) 

Here we have dropped the average pair momentum label and all primes on the momenta. 
Here also, Q = 2q (the relative momentum variable used in the analysis of pion correla- 
tions) and A is a fit parameter often called the chaoticity parameter. The matrix of radius 
parameters, [R 2 ], is the following real, symmetric matrix: 
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[R 2 



p2 d2 p2 \ 



p2 e>2 d2 



p2 p2 p2 



(17) 



in the Bertsch-Pratt parameterization. For pions, the Koonin-Pratt equation is a Fourier 
cosine transform that may be inverted analytically [[/J to give the source function directly in 
terms of the correlation function: 



S(f) = j^t- 3 J d 3 Q cos (Q ■ r) (C(Q) - l) . 



(2*) 

Inserting Eq. (|TB|) into this expression yields a Gaussian source function: 

A 



(18) 



S{r) = 

(2 v ^F) 3 y'det [R 2 [ 

Taking the r —> limit of this source 



exp (-KirjlR 2 }^ 



(19) 



A 



= (20) 
(2^)3^det [R 2 ] 

and inserting this into the equation for the average phase-space density, Eq. ( |14D or (|I5D 
yields the result in Ref. 0. Actually our result is more general than those in Refs. [0,|J as 
those results only apply either to diagonal [R 2 ] or [R 2 ] with R 2 ot ^ 0. One should note that 
the f = intercept of the source function has units of an effective volume. 



IV. PHASE-SPACE OCCUPANCY AND ENTROPY 

We can now estimate (/) from our calculation of the space-averaged phase-space density, 
(f(y,PT)) or (f(p)). In @, the authors argue that (/) can be estimated via 

/ d 3 p (f(p)f 

(/) = J -r- • (21) 

J d 3 p (f(p)) 

Given the current state of available correlation data, the uncertainty in the extracted sources 
can be greater than 50% of the extracted source value. Therefore, in this expression the 
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uncertainty in the result is dominated by the uncertainty in the underlying extracted source. 
Indeed, since the relative uncertainty in (f(p)) is nearly 5Sp(r — > 0)/Sp(r — > 0), we find 
S (f) I (f) ~ 3SSf(r — > 0)/Sp(r — > 0) for sources with a strong p dependence. Thus, one 
easily finds uncertainties greater than the values themselves. Clearly an alternative is needed 
that has a smaller dependency in the error on the source. 

Instead of (gTJ), we propose the following method for evaluating (/): 

1 r d 3 N 

M = 7^f*>VV»W (22) 

and for azimuthally symmetric systems, 

2tt f . r/-.V 



r d 2 N 

/ dydp T p T E (f(y,p T )) — — . (23) 



N part J 2ndydp T pT 

Using either of these expressions, the relative uncertainty goes like one factor of SSp(r 
0)/Sp(r — > 0). One can see (|22] ) (or (p3|)) by beginning with the definition of (/): 

/ d 3 rd 3 pf{f,p) 

(f) 



d 3 r d 3 p f(r,p) 

]d 3 p((f(p))Jd 3 rf(r,p)^) (24) 

N pai t/r 

which immediately gives Eq. (p2j). 

One might wonder if, through similar considerations, we may be able to improve on the 
calculation of the entropy per particle given in Ref. [0 or 0. In short, we do not believe so. 
Consider the entropy for a gas of fermions (top) or bosons (bottom): 

/ d 3 rd 3 p [/ln/±(lT/)m(lT/)] 

S/iV part = ~ l 7—— • (25) 

/ d 3 rd 3 pf 

To arrive at their expression for the entropy per particle, the authors of neglected the 
spatial dependence of f(r,p) and found: 
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S/N pavt = 

J d 3 p [(f(p)) In (f(p)} ± (1 =F (f(p))) ^ (1 T (f(p)))\ (26) 

d 3 p </(p)> ' 



The authors of || go one step farther and neglect the momentum variation of f(f, p) to obtain 
a result entirely dependent on the phase-space occupancy (/). Clearly the entropy should 
get larger contributions from regions of coordinate space where the phase-space density is 
small however both Eq. ( p6| ) and the analogous result from Ref. || do not reflect this. The 
source function's r dependence does give us information about the spatial dependence of the 
phase-space density however it is not clear how one might use this information to obtain a 
better estimate of the entropy. 

Now we present sample calculations of (/) and S/N paTt from two different experiments. 
The first calculation uses negative pion correlations measured from Pb+Pb collisions at 
158 GeVA from the CERN-SPS experiment NA49 HI3-15H- The second calculation uses 



proton correlations from the 14 N+ 27 A1 reaction at 75 MeVA measured at the Michigan 
State University NSCL |1(J. These calculations will show the general applicability of the 
source imaging and the superiority of Eq. ( p2|) over Eq. fl2~l~|). 

For the first calculation, we use the space-averaged tt~ phase-space densities of Ferenc et 
al. [0], extracted from experiment NA49 [p~3| p~5f| . In this calculation we could have used an 



extracted source, but we can also get the r' — > source intercept from the radius parameters 
of a Gaussian fit to the pion correlations, as we saw in section |111 Bj. Since this is exactly 



what is used in Ref. ||, we use their radii. We estimate the n spectrum as the product of 
the 7r~ rapidity and Pt distributions: 



0-9—3 — [TeMijn^ + Teff(y))] exp -- * — — (27) 



dydp T p T dy y T eS (y) 

Here, dN^/dy is the negative hadron rapidity distribution from Ref. and the factor 



of 0.9 accounts for the fraction of negative hadrons that are actually pions. In (|27|), the 
Pt distribution is parameterized by a rapidity dependent effective temperature T e s(y) and 



the actual values of this effective temperature are obtained from [ 15] . Using (p2|) , we find 

11 



(/) = 0.19 ± 0.06 while using (|T]), we find (/) = 0.14 ± 0.08. Both results are consistent, 
however the result from ( p2"D has a 25% smaller uncertainty. Using Eq. (pop, the entropy per 
pion can also be estimated as S/N n - = 3.9 ± 1.8. 



The numbers extracted from the sources from the NSCL pp data in [[ilj are more dra- 
matic. Since the proton spectrum for this reaction is only available at a few angles, we follow 
Ref. and approximate it with the thermal distribution: 

dN v 1 . , 

JL (x (28) 

dp z~ l exp (p 2 /2mT) + 1 

Here the normalization constant is determined by normalizing to the number of participants 
(9 protons), T is the fitted temperature of 10.2 MeV, and z is set from the requirement 
of maximum entropy giving z ~ 1.1. Since the source is given only at three fixed pair 
momenta, we follow Ref. and average the results to obtain the zero intercept S ave (r — > 
0) = 15.5 ± 2.5 x 10~ 4 fm~ 3 . Using (f2|), we find (/) = 0.27 ± 0.12 while using (0), we 
find (/) = 0.25 ± 0.04 - a factor of 3 improvement in the uncertainty. For the entropy 
per proton, we estimate S/N p = 2.7 ± 0.7. These results differ slightly from those in |7j 
simply because the authors of J7| set the error on the source to zero. Using the more 
accurately determined intercepts from Ref. ||, 5 ave (r — > 0) = 18.7 ± 1.1 x 10~ 4 fm~ 3 , we 
find (/) = 0.30 ± 0.02 from (gD and (/) = 0.30 ± 0.05 from (gl]). For the entropy we find 
S/Np = 2.45 ± 0.21. Note the substantial improvement from a better determination of the 
source intercept. Nevertheless, the result from Eq. (|22]) is still a factor of 2 improvement 
over (p|). 



V. CONCLUSIONS 

The phase-space density is an important, fundamental observable that can provide insight 
into the underlying dynamics of the nuclear reactions and it may be extracted from sources 
imaged in heavy-ion reactions. This extraction relies on an appropriate definition of the 
effective phase-space density at freeze-out, since the true freeze-out density is not a uniquely 
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defined concept. We have shown that the definition in Ref. |4[] is entirely consistent with the 
source function obtained by imaging J7J. We have also shown how the extraction of the space- 
averaged phase-space density from imaged sources is a generalization of Bertsch's result for 
pion correlations. Finally, we have provided a formula for the calculation of the average 
phase-space occupancy. This formula is less sensitive to the uncertainties of the source 
functions than others in the literature. We believe that the source imaging method will be 
useful for extracting space-averaged phase-space densities from data of future relativistic 
heavy-ion experiments. 
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